arXiv:1508.01525v2 [astro-ph.CO] 26 Feb 2016 


MNRAS 457, 3024-3035 (2016) 


1 March 2016 


Compiled using MNRAS IATeX style file v3.0 


Large-scale mass distribution in the Illustris simulation 

M. Haider^*, D. Steinhauser^, M. Vogelsberger^, S. Genel^f, V. Springe^’® 

P. Torrey^’® and L. Hernquist^ 

^ Institut fiir Astro- und Teilchenphysik, Universitdt Innsbruck, Technikerstrafie 25/8, A-6020 Innsbruck, Austria 
‘^Kavli Institute for Astrophysics and Space Researeh, Massaehusetts Institute of Teehnology, Cambridge, MA 02139, USA 
^Department of Astronomy, Columbia University, 550 West 120th Street, New York, NY 10027, USA 
"^Heidelberg Institute for Theoretieal Studies, Sehloss-Wolfsbrunnenweg 35, D-69118 Heidelberg, Germany 

^ Zentrum fur Astronomie der Universitdt Heidelberg, Astronomisches Recheninstitut, Monchhofstr. 12-14, D-69120 Heidelberg, Germany 
^ TAPIR, California Institute of Teehnology, Maileode 350-17, Pasadena, CA 91125, USA 
^Harvard-Smithsonian Center for Astrophysies, 60 Garden Street, Cambridge, MA 02138, USA 

Accepted 2016 January 8. Received 2015 December 21; in original form 2015 August 6 


ABSTRACT 

Observations at low redshifts thus far fail to account for all of the baryons expected 
in the Universe according to cosmological constraints. A large fraction of the baryons 
presumably resides in a thin and warm-hot medium between the galaxies, where they 
are difficult to observe due to their low densities and high temperatures. Cosmological 
simulations of structure formation can be used to verify this picture and provide 
quantitative predictions for the distribution of mass in different large-scale structure 
components. Here we study the distribution of baryons and dark matter at different 
epochs using data from the Illustris simulation. We identify regions of different dark 
matter density with the primary constituents of large-scale structure, allowing us to 
measure mass and volume of haloes, filaments and voids. At redshift zero, we find that 
49 % of the dark matter and 23 % of the baryons are within haloes more massive than 
the resolution limit of 2 x 10^ Mq. The filaments of the cosmic web host a further 
45 % of the dark matter and 46 % of the baryons. The remaining 31 % of the baryons 
reside in voids. The majority of these baryons have been transported there through 
active galactic nuclei feedback. We note that the feedback model of Illustris is too 
strong for heavy haloes, therefore it is likely that we are overestimating this amount. 
Categorizing the baryons according to their density and temperature, we find that 

17.8 % of them are in a condensed state, 21.6 % are present as cold, diffuse gas, and 

53.9 % are found in the state of a warm-hot intergalactic medium. 

Key words: galaxies: haloes - cosmology: dark matter - large-scale structure of 
Universe. 


1 INTRODUCTION 

The large-scale structure of the Universe is determined by 
the dominant dark matter, which makes up 26 % of to¬ 
day’s mass-energy content of the Universe. As dark matter 
is not detectable directly, the large-scale structure must be 
inferred from the baryons. Recent analysis (Planck Collab¬ 
oration 2015) of the cosmic microwave background fluctu¬ 
ations by the Planck mission found that baryons amount 
to — 4.8 % of the present-day critical density of the 
Universe. This corresponds to a mean baryon density of 
Pe = 4.1 X 10“^^gcm“^. An independent measurement of 
the baryon content of the Universe is possible through the 
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theory of big bang nucleosynthesis, as it allows one to in¬ 
fer the baryon content from the abundance ratio of light 
elements. The measurements of the abundance ratio of deu¬ 
terium to hydrogen (D/H) (Kirkman et al. 2003) are in very 
good agreement with the baryon content derived from the 
cosmic microwave background fluctuations. At higher red- 
shifts, the baryons were not fully ionized. This allows their 
detection in denser regions through Hi absorption lines in 
the spectrum of background quasars, the so-called Lyman 
alpha (Lya) forest. The amount of gas needed to produce 
the observed absorption lines in the Ly a forest is in good 
agreement with the cosmic baryon fraction (Weinberg 
et al. 1997). 

At low redshifts, however, observations fail to account 
for all the baryons detected at higher redshifts (Fukugita & 
Peebles 2004; Bregman 2007). While X-ray observations of 
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galaxy clusters do find baryon contents close to the primor¬ 
dial value (Simionescu et al. 2011), and Ovii absorption line 
measurements at X-ray energies hint at large reservoirs of 
hot gas around spiral galaxies (Gupta et ah 2012), approxi¬ 
mately 30 % of the baryons are still not accounted for in the 
local Universe (Shull et al. 2012). The baryons are not only 
missing on cosmological scales, but the baryon-to-dark mat¬ 
ter ratio in galaxies also falls short of the primordial ratio 
(Bell et al. 2003; McGaugh et al. 2010). 

Based on hydrodynamical simulations (Dave et al. 
2001), we expect that a significant amount of the baryons 
are hidden in the state of a thin warm-hot intergalactic 
medium (WHIM) with temperatures between 10^ and 10^ 
K. Through Hi and O VI surveys (Danforth & Shull 2008), 
it is possible to probe the colder parts of the low-redshift 
intergalactic medium. However, the gas is too diffuse and 
the temperatures are insufficient to be detectable with the 
current generation of X-ray satellites (Kaastra et al. 2013), 
except in the densest regions between galaxy clusters (Nicas- 
tro et al. 2005, 2013). The question of how the baryons are 
distributed and whether they follow the filaments of the cos¬ 
mic web is not well constrained observationally, and con¬ 
sequently, cosmological simulations of structure formation 
are important tools for making predictions about the large- 
scale distribution of baryons (Gen & Ostriker 1999, 2006; 
Dave et al. 1999, 2001; Smith et al. 2010). In this paper, 
we investigate the distribution of baryons and dark matter 
using data from the Illustris simulation (Vogelsberger et al. 
2014b). In Section 2, we give a short overview of the sim¬ 
ulation and discuss the methods we used. Our results are 
presented in Section 3. Specifically, in 3.1, we investigate 
the amount of matter inside haloes and compute the baryon 
fraction for haloes of different masses. In Section 3.2, we ex¬ 
amine the distribution of matter with respect to the dark 
matter density. This allows us to decompose the simulation 
volume into haloes, filaments and voids, and measure the 
mass, volume and redshift evolution of these components. 
In Section 3.3, we look at the distribution of the baryons in 
the temperature-density space and compute the mass frac¬ 
tion and redshift evolution of the condensed, diffuse, hot and 
WHIM phases. We discuss the results in Section 4 and close 
with a summary in Section 5. 


2 SIMULATION METHODS 

We analyse data from the Illustris simulation, which is a 
cosmological hydrodynamics simulation of galaxy formation 
in a (106.5 Mpc)^ volume [corresponding to a (75 Mpc/h)^ 
periodic box]. The gas and dark matter were evolved using 
the AREPO moving-mesh code (Springel 2010). In addition 
to gravity and hydrodynamics, the simulation includes sub¬ 
grid models for star formation, active galactic nuclei (AGN) 
feedback and gas cooling [see Vogelsberger et al. (2013) and 
Torrey et al. (2014) for details on the implemented feed¬ 
back models and the parameter selection]. Dark matter and 
gas are represented by 1820^ resolution elements each, re¬ 
sulting in an initial mass resolution of 1.26x10® M© for 
gas and 6.26x10® M© for dark matter. The simulation uses 
a standard AGDM model with H = 70.4, Qo = 0.2726, 
= 0.7274 and 0.0456. It was started at redshift 

z = 127 and evolved all the way to the present epoch. The 


Illustris simulation produced a galaxy mass function and a 
star formation history which are in reasonably good agree¬ 
ment with observations. Remarkably, it also yielded a re¬ 
alistic morphological mix of galaxies. An overview of basic 
results is given in Vogelsberger et al. (2014b), Vogelsberger 
et al. (2014a) and Genel et al. (2014), and recently, all of 
the simulation data has been made publicly available (Nel¬ 
son et al. 2015). 

In Section 3.1, we use the halo catalogue of the Illus¬ 
tris simulation, which has been generated using the SUB¬ 
FIND halo-finding algorithm (Springel et al. 2001). SUBFIND 
searches for gravitationally bound overdensities and com¬ 
putes the dark matter, gas and stellar mass which is bound 
to a halo or subhalo. In the subsequent analysis, we com¬ 
pute the density fields of gas, stars, baryonic matter and 
dark matter on a uniform grid with 1024^ cells. To this end, 
the dark matter and gas particles are mapped conservatively 
to the grid using a smoothed particle hydrodynamics (SPH) 
kernel technique. For the dark matter, we set the smoothing 
length to be the radius of a sphere containing the 64 near¬ 
est dark matter particles, as is done also in SUBFIND. The 
gas cells are mapped using three times rc = (3V[.gii/47r)^/®, 
where is the volume of a Voronoi cell. The stellar and 
black hole particles have been mapped using the particle in 
cell method with a nearest grid point assignment, as they 
are not tracers of an extended density field. We should note 
that the densities given in this analysis are thus average den¬ 
sities over the volume of individual cells, which are 104 kpc 
on a side. 

The baryon densities and temperatures used in Sec¬ 
tion 3.3 have been computed directly from the Voronoi cells 
of the simulation and do not involve any additional smooth¬ 
ing. For comparison, we also use a variant of the Illustris 
simulation where cooling, star formation and feedback have 
been switched off. This simulation, which we will refer to 
as non-radiative run, has only half the resolution of the full 
physics simulation. 


3 RESULTS 

In Section 3.1, we investigate the mass fraction of baryons 
and dark matter in haloes and compute the baryon-to-halo 
mass ratio. In Section 3.2, we examine the mass distribution 
with respect to the dark matter density. With the help of 
cuts in dark matter density, we divide the simulation volume 
into haloes, filaments and voids. In Section 3.3, we analyse 
the baryon distribution in temperature-density space. 

3.1 Mass in haloes 

In Table 1, we give the mass fraction of dark matter and 
baryons residing inside the identified SUBFIND haloes. The 
mass fractions are calculated using the mass of all parti¬ 
cles which are gravitationally bound to haloes. Only main 
haloes were considered, with the mass of genuine subhaloes 
being added to their parent halo. The table also shows the 
decomposition of the baryonic mass into gas and stars. We 
should note that the mass found in haloes depends on the 
resolution of a simulation, therefore Table 1 only gives the 
fraction of mass in haloes which can be resolved in the Il¬ 
lustris simulation. The resolution limit in terms of mass is 
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primarily given by the mass of a single dark matter parti¬ 
cle which is 6.26 x 10® M©. To detect gravitationally bound 
haloes, a minimum of about 30 particles are required. Thus, 
the threshold for halo detection in Illustris is approximately 
2 X 10® M 0 . A resolution study using a run with one-eighth 
of the mass resolution leads to similar results as the ones 
presented in Table 1. This suggests that the results are nu¬ 
merically robust. How much of the dark matter is expected 
to reside in haloes below the resolution limit also depends 
on the physical nature of the dark matter particle. Angulo & 
White (2010) find that for a 100 GeV neutralino dark mat¬ 
ter particle only 5-10 % of the dark matter should not be 
part of a clump at redshift zero, and that the most abundant 
haloes have masses around 10 ~® M©. 

We find that about half of the total dark matter mass 
in the simulation volume is contained within haloes, while 
only 21 % of the baryonic mass resides within haloes. This 
evidently means that in the Illustris simulation, the mass 
ratio between baryons and dark matter is on average lower 
in haloes than the primordial mixture, f^b/^dm- latter 
would be the expected value if baryons traced dark matter 
perfectly. The second row of Table 1 shows that nearly all 
the halo baryons are in haloes with a total mass higher than 
5x10^ Mq. Less massive haloes in Illustris consist mainly 
of dark matter and are hence truly dark. We should note 
though that a 10^ M© halo is made up of only 160 dark 
matter particles, and the results will be less accurate at this 
scale than for haloes with higher numbers of particles (Vo- 
gelsberger et al. (2014a), see Trent! et al. (2010) for a dis¬ 
cussion of the influence of particle number on the quality of 
halo properties). For haloes more massive than 5x10^ M©, 
we again see that the baryon/dark matter ratio is only half 
of the primordial ratio, and in the mass range between 10 ^^ 
and 10^® M©, it drops to about one-third. Towards the most 
massive haloes of cluster size, this rises again to about 50 % 
and beyond of the expected primordial value. 

Comparing these mass fractions bound in haloes with 
the results of the non-radiative run of the Illustris initial con¬ 
ditions, we find that the amount of dark matter in haloes is 
very similar. However, in the non-radiative run, the fraction 
of baryons in haloes is much higher and reaches about 40 %, 
which is roughly double the value of the full physics run. As 
the overall dark matter structure is very similar for the full 
physics and the non-radiative run, the difference in baryonic 
mass is a consequence of the feedback processes included in 
the full physics simulation (see also Sections 3.2 and 4). 

The mass fraction in haloes also depends on the defini¬ 
tion of halo mass. The values presented in Table 1 are the 
total gravitationally bound masses. Another common way 
to define the mass of haloes is to count the mass inside a 
sphere whose radius is set so that the mean enclosed density 
equals a reference density. If we measure the mass around 
the primary haloes of the SUBFIND friend-of-friend groups, 
and use as a reference density 200 x pmean (where pmean 
denotes the mean matter density of the Universe) we find 
a total mass fraction of 48.1 % in haloes. Using the higher 
reference density 500 x this value reduces to 26 %. We 
also want to note that for high-mass haloes, it makes a dif¬ 
ference if the mass includes the contribution of subhaloes or 
not. In Illustris, 19.2 % of the haloes are subhaloes. They 
make up for 13.7 % of the total mass and 29.7 % of the 
stellar mass. Thus, if we neglected the contribution of sub- 



Figure 1. Baryon to total matter ratio inside central haloes of 
the Illustris simulation divided by the primordial ratio 
(solid black line). The x-axis gives the total mass of haloes inside 
a sphere with a mean density of 500 times the critical density. The 
dashed black line shows the contribution of stars, the red line the 
contribution of gas warmer than 10^ K, the blue line gas colder 
than 10^ K and the green line the amount of neutral hydrogen. 
The black squares and grey triangles represent observational data, 
taken from table 2 of McGaugh et al. (2010). 



Figure 2. Redshift evolution of the mass inside haloes. The val¬ 
ues for dark matter and baryons are normalized to the total dark 
matter and baryon mass in the simulation volume, respectively. 
The dashed lines show the contribution of gas and stars to the 
baryon fraction. 

haloes, the haloes with a mass higher than lO^^M© would 
be found to host only 4.4 % of the total matter, 4.9 % of the 
dark matter and 1.9 % of the baryons. 

We further illustrate the baryon-to-dark matter ratio 
in Fig. 1, where we plot the ratio of the baryon content 
Mbaryon/^total primordial baryon fraction 

against the halo mass. The baryon fraction is broken down 
into components of cold gas, hot gas, neutral hydrogen and 
stars, and the plot includes observational data compiled by 
McGaugh et al. (2010) for comparison. Only central haloes 
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Table 1. Mass fraction of dark matter and baryons inside haloes at 2 : = 0. The mass fractions are given with respect to the total mass 
inside the simulation volume. The first row shows the fractions for all haloes in the Illustris full physics simulation. The second row gives 
the mass contribution of haloes with a total mass higher than 5 X 10® Mq and the third and fourth rows show the respective values 
for haloes more massive than 10^^ and 10^^ Mq, respectively. The last row gives the mass fraction for all haloes of the non-radiative 
Illustris run (no star formation, feedback or cooling). 



% of total 
dark matter mass 

% of total baryonic 
baryons gas 

mass 

stars 

% of total 

mass 

all haloes 

50.4 % 

21.3 % 

14.7 % 

6.6 % 

45.5 % 

Mtot > 5 x109Mo 

44.6 % 

21.3 % 

14.7 % 

6.6 % 

40.7 % 

Mtot > 

29.3 % 

9.8 % 

4.8 % 

5.0 % 

26.0 % 

Mtot > 1013 m© 

20.1 % 

7.2 % 

4.0 % 

3.2 % 

18.0 % 

Mtot>10i^M© 

11.8 % 

6.0 % 

3.7 % 

2.3 % 

10.8 % 

all haloes, non-radiative run 

49.8 % 

39.1 % 

- 

48.0 % 


(those which are not subhaloes) have been used for this plot 
and only the mass inside i?500c (^^e radius inside of which 
the mean density is 500 times the critical density) is taken 
into account. A resolution study using data from a run with 
one-eighth of the mass resolution produced results very sim¬ 
ilar to Fig. 1. 

In the regime of galaxy clusters. X-ray observations 
find a baryon to dark matter ratio between 70 and 100 % 
of the primordial value (Vikhlinin et al. 2006). For clus¬ 
ter sized haloes, Illustris finds a ratio around 40-50 %, but 
only small clusters are present in the limited simulations 
volume. This disagreement with observations has already 
been noted (see Fig. 10 in Genel et al. 2014) and has been 
attributed to a too violent radio-mode AGN feedback model. 
The radio-mode feedback model is used when the black hole 
accretion rates are low. It is most effective at removing gas 
from the haloes in the mass range from 10^^'^ to 10 ^^'^Mq, 
as at lower halo masses, the black holes are less massive 
and at higher halo masses, the potential wells are deeper. 
The model and the used parameters result in too large mass 
outflows from group-sized haloes and poor clusters (see Sec¬ 
tion 4 for a more detailed discussion of the model). As a 
considerable amount of mass is within group-sized haloes 
and clusters, the results about the distribution of baryons, 
which we present in the following sections, will be affected 
by the mass outflows. 

Between 10^° Mq and 10^^ Mq, the baryon fraction is 
close to the primordial value In this mass regime, 

the radio-mode AGN feedback plays only minor role as the 
black holes are less massive. The applied supernova feedback 
model is removing gas from the star-forming phase through 
galactic winds. However, the velocity of these winds is lower 
than the escape velocity and thus the supernova feedback 
model does not remove gas from the haloes. We find that 
most of the baryons in this mass range are at relatively low 
temperatures, and a large fraction of the gas is in neutral 
hydrogen. We should note though, that Illustris does not 
simulate the effects of the stellar radiation field on the gas, 
and therefore we cannot predict the ionization state exactly. 
The baryon fraction we find is higher than the fraction of 
detected baryons in McGaugh et al. (2010). However, these 
observations should be thought of as lower limits. A large 
fraction of the baryonic content of galaxies is expected to re¬ 
side within the circumgalactic medium, which is subject to 
large observational uncertainties and is depending on model 
assumptions. Recent studies (Tumlinson et al. 2013; Bor- 
doloi et al. 2014; Werk et al. 2014) suggest that galaxies can 


have substantial amounts of cold gas in their circum-galactic 
medium. In particular, we note that the baryon fraction we 
find in Illustris is compatible with the allowed range deter¬ 
mined by Werk et al. (2014) (see their Fig. 11). The high 
baryon fraction in this mass range might also be connected 
to the finding that the galaxy stellar mass function in Il¬ 
lustris is too high for low mass galaxies (Genel et al. 2014; 
Vogelsberger et al. 2014a). We also recall that arepo leads 
to more cooling than SPH-based codes for the same physics 
implementation (Vogelsberger et al. 2012). 

In Fig. 2, we show the redshift evolution of the fraction 
of total mass, dark matter and baryon mass bound inside 
haloes. We see that the total mass inside haloes is mono- 
tonically increasing with time, reaching 45.5 % at redshift 
z = 0. Between redshift z — 2 and 1, the mass gain for the 
baryonic halo component is less pronounced than for dark 
matter, and from redshift z = 1 to 0 the fraction of baryons 
inside haloes even decreases. The onset of this decrease coin¬ 
cides with the radio mode AGN feedback becoming impor¬ 
tant (see Sijacki et al. 2015, for the evolution of the black 
hole accretion rate in the Illustris simulation). The figure 
also shows that the dominant contribution to the baryonic 
halo budget is from gas at all redshift. While the fractional 
contribution of stars to the total baryon budget increases 
with time, it reaches only a final value of 6.6 % at redshift 


3.2 Mass distribution in different dark matter 
density environments 

The formation of large-scale structure is dominated by dark 
matter. Therefore, in this section we analyse the mass dis¬ 
tribution with respect to different dark matter density en¬ 
vironments. We compute the dark matter density through 
mapping the Illustris data on to a (1024)^ grid covering the 
whole simulation volume. In Fig. 3 (a) we show a plot of the 
dark matter density and in Fig. 3 (b) of the baryon density 
in a slice through the simulation at redshift zero. The slice 
has an extent corresponding to the full 106.5 x 106.5 Mpc, 
and a thickness of one cell (104 kpc). For both the dark 
matter and the baryon slice we use a logarithmic colourmap 
with a range of 4.5 dex. The colourmap boundaries for the 
baryons have been set an order of magnitude lower than for 
the dark matter. 

We note that the baryons appear more extended than 
the dark matter distribution, especially around high dark 
matter density regions. It is interesting to compare this to 
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Fig. 4, where we show the corresponding baryon density of 
the non-radiative run. Clearly, in the absence of feedback, 
the baryons are tracing the dark matter well on large scales. 

Comparing the baryon distributions at different red- 
shifts (Figs. 5 and 3), one can see that the extended baryon 
distributions around high dark matter density areas are less 
prominent at 2 : = 0.5 and virtually absent at 2 : = 2. The 
diameter of these regions is also increasing from z = 0.5 
to 0. We conclude that the extended gas shells originate in 
gas which has been expelled from their host haloes due to 
radio-mode feedback, which becomes active at low redshift. 
A rough estimate from the diameter of these shells at differ¬ 
ent redshifts suggests that their expansion speed is between 
500 and 1000 comoving kms“^. 

3.2.1 Mass distribution according to dark matter density 

In Fig. 6, we present the fraction of the baryons at a given 
baryon and dark matter density. We binned the data in 
baryon and dark matter density using 200 x 200 bins. 
The diagonal line gives which is where the 

baryons would lie if they traced dark matter perfectly. The 
colourmap shows the contribution of a bin to the total 
baryon mass. We find that there is a significant fraction of 
mass at high dark matter and high baryon densities, which 
corresponds to gas in haloes. A large fraction of baryons is 
at intermediate dark matter densities around the primordial 
dark matter density f^dmPcrit- However, there are also many 
baryons at the lowest dark matter densities, the region we 
labelled ‘voids’. 

It is interesting to observe the differences between the 
full physics and the non-radiative runs in Fig. 7, which shows 
the distribution of mass according to the dark matter den¬ 
sity. We binned the data in dark matter density and then 
measured the mass contribution of a bin to the total bary- 
onic or dark matter mass in the simulation volume. We see 
that the dark matter distribution in (b) is very similar for 
the full physics and the non-radiative run. However, there 
are notable differences for the baryons. For the full physics 
simulation, there is less mass at higher dark matter densities 
and significantly more mass at the lowest dark matter den¬ 
sities. This again is due to the feedback processes present 
only in the full physics simulation, which expel matter from 
the haloes out into dark matter voids. 

3.2.2 Dark matter density and the cosmic web 

The primary constituents of the large-scale structure differ 
in their dark matter density. Therefore we can use the dark 
matter density itself as a rough proxy to split up the simula¬ 
tion volume into regions corresponding to haloes, filaments 
and voids. We thus define the categories haloes, hlaments 
and voids by assigning to them a certain dark matter den¬ 
sity range^. By summing the grid cells which fall into the 
respective dark matter density ranges we can measure the 
mass and volume of those regions. We show the resulting 
mass and volume fraction and the density range used for 
this classification in Table 2. The spatial distribution of the 

^ The dark matter densities we use in this analysis are average 
densities over a volume of (104 kpc)^. 


categories in Table 2 is shown in Fig. 8. We show the same 
slice as in Fig. 3, and use a white colour for those cells that 
belong to the corresponding category. 

The ‘haloes’ are defined to have a dark matter density 
higher than I5pcrit- This density threshold has been cho¬ 
sen so that the halo region contains approximately the same 
amount of dark matter as the haloes of the halo finder. Thus, 
by construction, we find 49 % of the dark matter in haloes. 
But also the amount of baryons, which corresponds to 23 %, 
is in good agreement with Table I. The volume fraction of 
the haloes is 0.16 %. The chosen threshold is also consistent 
with the resolution of the grid on which we calculate the den¬ 
sity and the mass resolution of Illustris: a halo of 2 x 10^ M© 
would result in an overdensity of I5xp(,j.j^, if all of its mass 
falls into a single grid cell. According to the SUBFIND halo 
catalogue, haloes smaller than 2x10^ M© still host 4.8 % of 
the dark matter. However, many of those haloes will be sub¬ 
haloes and thus contribute their mass to the halo category. 
If we define the halo category to have densities higher than 
30pcrit? T hosts 41.7 % of the dark matter and 

21.6 % of the baryons, and all of the stars. This density cut 
at 30pcrit corresponds to the mass of a 4.2 xIO® M© halo 
in one cell. If we compare this value with the second row in 
Table I > 5 x 10^ M©) we again find that the value 

obtained using the dark matter density method is consistent 
with the value of the SUBFIND halo catalogue. We want to 
emphasize that the threshold density for the haloes is de¬ 
pending on the resolution of the grid, as it is the average 
density in one grid cell, and not the physical density at the 
outskirts of haloes. 

The ‘hlaments’ region, which we dehned to have dark 
matter densities between 0.06 to 15 Pcrit’ hosts 44.5 % of the 
dark matter and 46.4 % of the baryons. Its volume makes up 
21.6 % of the total simulation volume. The hlaments span a 
large density range in our dehnition. At the higher density 
end of this dehnition, between 5 and 15 Pcrit? ^he mass is 
mainly in ring-like regions around the bigger haloes, and no 
hlamentary structure is visible in this density regime. This 
circumhalo region hosts about 11.1 % of the dark matter and 
3.3 % of the baryons. The density boundary to the voids is 
motivated by Fig. 6, as the large bulk of mass in the denser 
parts of the hlaments extends down to this value. Also, in¬ 
specting Fig. 8, this density range seems to correspond to 
what one would visually label as hlaments. However, a dif¬ 
ferent density range could as well be justihed. 

The ‘voids’ (dark matter densities between 0 and 
0.06 Pcrit) contain only 6.4 % of the dark matter but 30.4 % 
of the baryons. The volume of the voids makes up 78.2 % of 
the total simulation volume. In that respect it is worthwhile 
to examine Figs. 10 and II, where we show the temperature 
and the metallicity of the gas. If we compare the temper¬ 
ature and metallicity maps with Fig. 3 and Fig. 8, we see 
that the majority of the baryons in voids is composed of 
warm to hot gas enriched by metals. Those baryons have 
most likely been ejected from the haloes through feedback. 
Because the ejected material has a higher temperature than 
the other baryons in voids, we can use the temperature to 
discriminate between baryons naturally residing in voids and 
baryons which have been transported there. We define an ad¬ 
ditional ‘ejected material’ region in Table 2 which is defined 
as having a temperature higher than 6 x 10^ K in addition 
to the dark matter density cut. With 23.6 % of the total 
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(a) dark matter (b) baryons 


Figure 3. Dark matter and baryon density in a thin slice at 2 : = 0. The slice covers the whole (106.5 Mpc)^ extent of the simulation 
and has a thickness of 104 kpc (1 cell). 


Table 2. Dark matter mass, baryonic mass and volume fraction in haloes, filaments and voids at z=0. The categories have been defined 
through dark matter density ranges. We also added a category ‘ejected material’ which corresponds to baryons inside the ‘voids’ region 
which have a temperature T> 6 x 10^ K. The spatial regions to which these dark matter density regions correspond to are shown in 
Fig. 8. 


dark matter density % of total % of total % of total % of total 


component region (p^j-it) dark matter mass baryonic mass mass volume 


haloes 

> 15 

49.2 % 

23.2 % 

44.9 % 

0.16 % 

filaments 

0.06 - 15 

44.5 % 

46.4 % 

44.8 % 

21.6 % 

voids 

0 - 0.06 

6.4 % 

30.4 % 

10.4 % 

78.2 % 

ejected material 
inside voids 

0 - 0.06 

2.6 % 

23.6 % 

6.1 % 

30.4 % 


baryons, this ‘ejected material’ region is responsible for most 
of the baryons in dark matter voids. In Fig. 9, the spatial 
region corresponding to the ejected material is plotted; note 
that it fills about 40 % of the voids. We should note though, 
that the ejected mass most likely heats some of the baryons 
already present in the voids. Therefore, we have probably 
overestimated the ejected mass in voids. However, through 
following the redshift evolution of the mass in voids we can 
give an estimate of the associate uncertainty, as we discuss 
below. We note that our findings for the volume fractions 
are generally in good agreement with simulations by Can- 
tun et al. (2014). 


3.2.3 Redshift evolution of matter and metals in haloes, 
filaments and voids 

By applying the same dark matter density cuts at different 
redshifts, we can study the time evolution of the values re¬ 
ported in Table 2. This is done in Fig. 12, where we show 
how the baryons and dark matter divide into haloes, fila¬ 
ments and voids as a function of time. In Fig. 12 (a) we see 
that, starting at redshift z = 2, feedback begins to efficiently 
remove gas from haloes. At first, this only slows down halo 


growth, but after a redshift of z = 1, it reduces the amount 
of baryons in haloes. In Fig. 12 (b) we see that the dark mat¬ 
ter haloes, unaffected by feedback, continue to grow at the 
expense of the filaments. At high redshifts, the dark matter 
was distributed homogeneously with a density of ^dmPcrit? 
and thus falls into the ‘filament’ category. The underdense 
regions of the voids were only created as matter from less 
dense regions was pulled into denser regions. Thus the frac¬ 
tion of dark matter in voids is increasing from z = 6 to 2. 
After z = 2, the amount of dark matter in voids is slowly de¬ 
creasing due to accretion on to filaments. The baryons show 
a similar behaviour from z = 6 to 2. However, starting at 
a redshift of z = 2 ejected material is also transported into 
the voids, thus increasing their baryonic content. We see in 
Fig. 12 (a) that the mass increase of the ‘ejected material’ 
is higher than the mass increase of the ‘voids’. The most 
likely explanation is that the ejected mass heats gas already 
present in the voids, which means that we overestimate the 
mass of the ejected material with our density and tempera¬ 
ture cut. If we assume that in the absence of feedback the 
baryons would show the same relative decrease from z = 2 
to 0 as the dark matter, we would need to correct the value 
of the ejected material down to 20 %. 
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Figure 4. Baryon density of the non-radiative simulation (no 
star formation, feedback or cooling) at z = 0. The same slice as 
in Fig. 3 is displayed. 


Illustris can also be used to probe the metal content in 
gas and stars. Star particles are stochastically formed when 
gas reaches densities above a threshold value A star 
particle is modelled as a single-age stellar population, for 
which the mass and metal return will be computed at every 
time-step. This material is then distributed over neighbour¬ 
ing gas cells (see Section 2.2 and 2.3 in Vogelsberger et al. 
2013). In Fig. 13, we show the evolution of the metals in 
gas and stars normalized to the total metal mass at z = 0. 
Additionally, the figure shows which fraction of the metals 
in gas is residing in haloes, filaments or voids. We find that 
at z = 0, 36 % of the metals are locked up in stars and 64 % 
of the metals are in gas. Considering only the gaseous com¬ 
ponent, half the metals are within haloes and 28 % reside in 
the filaments. The remaining 22 % are located in voids. The 
average metallicity of the stars is 1.49 solar metallicities at 
z — 0, while the halo gas has about 0.37 solar metallicities 
on average. The average metallicity in filaments and voids 
is roughly 0.1 times the solar value. 


3.3 Baryonic temperature—density relation 

An alternative way to look at the distribution of baryons is 
to analyse them according to their density and temperature. 
This is more directly relevant to observations, as density and 
temperature are the important variables for emission and 
absorption mechanisms. In Fig. 14, we show the contribu¬ 
tion of gas to the total baryonic mass in a temperature versus 
baryon density histogram. For this plot we directly use the 
Voronoi cell densities of the gas instead of the averaged den¬ 
sities we used in Section 3.2. We divide the baryons accord¬ 
ing to the same classification as has been used in Dave et al. 
(2001): diffuse gas having p < lOOOpcrit^b ^ < 10^ K, 
condensed gas with p > lOOOpcrit^b ^ ^ warm- 

hot gas with temperatures in the range 10^ <T < WK 


and hot gas with temperatures above 10^ K. We refer to the 
warm-hot gas also as WHIM. 

We find that 21.6 % of the baryons are in the form of 
diffuse gas, located mainly in the intergalactic medium. The 
tight relation between temperature and density is due to the 
interplay of cooling through adiabatic expansion and pho¬ 
toionization heating (see also the discussion in Section 3.3 
of Vogelsberger et ah, 2012). The condensed gas amounts to 
11.2 %. Most of the condensed gas is in a horizontal stripe 
around 10^ K. As photoionization heating becomes less dom¬ 
inant but cooling time-scales shorten at higher densities, gas 
cools effectively down to this temperature. Since there is no 
metal and molecular line cooling, 10^ K represents an effec¬ 
tive cooling floor for the gas in this phase. The upward ris¬ 
ing slope which extends from the ‘condensed’ region into the 
‘WHIM’ region is due to an effective equation of state used 
for gas exceeding the density threshold for star formation 
(Springel V Hernquist 2003). The warm-hot medium makes 
up for 53.9 % of the baryons, while 6.5 % of the baryons are 
in hot gas. 

The mass in these two categories corresponds to warm- 
hot shock-heated gas in haloes and filaments and material 
which has been ejected due to feedback from haloes (in¬ 
spection of Fig. 10 shows that the ejected material is in the 
temperature range defining the WHIM region). If all the 
23.6 % of the ejected material (see Table 2) had remained 
in haloes, this would have changed the warm-hot mass frac¬ 
tion to 30 % and increased the condensed fraction to 34.8 % 
(plus 6.6 % in stars). 

The redshift evolution of the gas phases is given in 
Fig. 15, where we see that at high redshift, most of the gas 
has been in the form of a rather cold and diffuse medium. 
Starting at a redshift of z = 4, and more pronounced after 
redshift z — 2^ the WHIM phase is gaining more and more 
mass, and by redshift zero, ends up containing most of the 
baryons. 


4 DISCUSSION 

Comparing the values of Table 1 with Table 2, we see a 
very good agreement at z = 0 between the SUBFIND halo 
catalogue and the haloes defined by a dark matter density 
cut. This suggests that our method of measuring the mass 
using the average dark matter density in a cell of our grid 
works reasonably well. However, we see deviations of 5-7 % 
for the mass in haloes at redshifts higher than 1. The reason 
for these deviations is that at higher redshifts the haloes are 
less massive, and thus some fall below the resolution limit 
of our grid. 

Using the temperature-baryon density classification in 
Section 3.3, we find that 53.9 % of the baryons reside in 
the WHIM region. This is higher than the 30-40 % found 
in the work of Dave et al. (2001), or the work of Cen V 
Ostriker (2006), who reported between 40 and 50 % in the 
WHIM phase. This discrepancy is most likely due to the use 
of different feedback models. The importance of the feedback 
model is further underlined by the differences between the 
full physics and the non-radiative runs, which produce a 
nearly identical dark matter distribution but very different 
baryon distributions (see Fig. 7). In the full physics run only 
half the baryons are within haloes compared to the non- 
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Figure 5. Redshift evolution of the baryon density. The same slice as in Fig. 3 has been used. 
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Figure 6. Histogram of the contribution of different dark mat¬ 
ter and baryon density regions to the total baryonic mass. The 
dashed vertical lines indicate the boundaries of the dark matter 
density regions we defined in Section 3.2.1 to measure the mass 
in filaments, voids and haloes. The diagonal line indicates where 
a primordial mix of baryons and dark matter would lie. 


radiative simulation (see Table 1). This difference is even 
more significant when we consider that the non-radiative 
run contains no radiative cooling and therefore the baryons 
cannot condense as easily into haloes as in the full physics 
run. 

As stated in Section 3.2.1, it is interesting to note that 
we find 30 % of the baryons residing in voids. The reason for 
this is the strong radio-mode AGN feedback model, which is 
expelling baryonic matter from haloes and is responsible for 
nearly 80 % of the baryons ending up in voids. In a recent 
survey, Shull et al. (2012) find that 30 % of the baryons are 
missing in the low-redshift Universe. It would be tempting 
to explain these missing baryons with the baryons we find in 
voids. However, we have to be cautious in reaching this con¬ 
clusion, as it is not clear whether the radio-mode feedback 
model implemented in Illustris is reliable. An indication that 
the feedback model of Illustris is too strong comes from the 
gas content of massive haloes. Observations of galaxy clus¬ 


ter sized haloes find baryon fractions between 70 and 100 % 
of the primordial value (Vikhlinin et al. 2006). In contrast, 
we find (see Table 1) that the Illustris simulation produces 
a baryon fraction of only 50 % for haloes more massive than 
lO^^M©. The low gas content of massive haloes in the sim¬ 
ulation has been attributed to the strong AGN radio-mode 
feedback (Genel et al. 2014). The radio mode is used when 
the accretion rate on to a black hole is low, and thus is 
strongest at low redshifts. Therefore, our results at low red- 
shifts would most likely change if a weaker AGN feedback 
model was used. 

Do the missing baryons in massive haloes explain the 
baryons we find in voids? The 149 haloes with a total mass 
higher than IO^^Mq contain 20.1 % of the total dark matter 
but only 7.2 % of the baryonic mass. Thus, 13 % of the 
total baryons are missing in these haloes. It therefore seems 
that the missing mass cannot explain the full 24 % of ejected 
material in voids. However, it is possible that material which 
is pushed out from the haloes might remove gas from the 
filaments. We should also note that the 24 % we quote is 
only the ejected mass inside voids. The actually ejected mass 
might be higher as some of the ejected material would still 
be in the filaments. 

Other simulations of structure formation which include 
AGN feedback have produced baryon fractions in massive 
haloes in agreement with observations [McGarthy et al. 
(2011), Le Brun et al. (2014), Schaye et al. (2014), Planelles 
et al. (2013) or Khandai et al. (2015)]. It is beyond the scope 
of this paper to give a detailed comparison of different ap¬ 
proaches. However, as the AGN feedback has a large im¬ 
pact on our results, we want to highlight some differences in 
the models. Most current cosmological simulations, includ¬ 
ing Illustris, include AGN feedback through modifications 
and extensions of the approach introduced in Springel et al. 
(2005). There, black holes are seeded when the hosting halo 
passes a threshold mass (5 x 10^° h~^MQ for Illustris), and 
then Eddington-limited Bondi-Hoyle accretion is assumed. 
As the actual accretion region cannot be resolved, and the 
medium there is usually a mixture of hot gas and cold clouds, 
the Bondi-Hoyle accretion is sometimes enhanced by a fac¬ 
tor a, which is invoked to compensate for the lack of reso¬ 
lution. In Illustris, this accretion enhancement parameter is 
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(a) (b) 

Figure 7. Baryon and dark matter mass distribution with respect to the dark matter density at 2 : = 0. In (a), it is shown how much 
baryonic mass resides in an environment with a given dark matter density. The y-axis shows the fraction of baryonic mass in a dark 
matter density bin to the total baryonic mass in the simulation volume. The solid line represents the data from the full physics run and 
the dotted line shows the values for the non-radiative run. The dashed line gives the contribution of stars to the baryons in the full 
physics run. In (b), we show the analogous plot for the dark matter mass. 



(a) voids (b) filaments (c) haloes 


Figure 8 . Spatial regions corresponding to the dark matter density ranges defined in Table 2. If a cell’s dark matter density is inside 
the corresponding region it is drawn in white, otherwise it is black. The same slice as in Fig. 3 has been used. 


set to a = 100. It is further assumed that a fraction Cr (set to 
0.2 in Illustris) of the rest mass energy of the accreted mass 
is converted into energy, and a fraction of this energy in turn 
couples to the surrounding material. In Illustris, a modifi¬ 
cation of the model by Sijacki et al. (2007) is used, where 
a distinction between a quasar mode and a radio mode has 
been introduced. The quasar mode is assumed to be active 
if the accretion rate is above 5 % of the Eddington rate, and 
the radio mode, is used below this threshold. In the quasar 
mode, a fraction Cf = 0.05 of the AGN luminosity is coupled 
thermally to cells around the black hole. In the radio mode, 
the AGN is assumed to inject energy into bubbles, whose ra¬ 
dius and distance to the black hole are computed according 
to Sijacki et al. (2007). The radio mode is assumed to have 
a higher coupling efficiency of em = 0.35. The radio mode 
bubbles are not formed continuously, but only after the mass 
of the black hole increased by a specified value (15 % in 
Illustris). If the fractional mass increase AM/M surpasses 


energy emCr AMc^ is put into a bubble. Thus, 
controls the frequency and energy of bubbles. 

Among the recent hydrodynamical simulations of 
galaxy formation mentioned above, Illustris is the only sim¬ 
ulation which makes a distinction between radio and quasar 
mode in this way. Planelles et al. (2013) also use different 
efficiencies for radio and quasar mode, but the energy is 
coupled to surrounding particles in the same way for both 
modes. The model of Booth & Schaye (2009), which is used 
in McGarthy et al. (2010), Le Brun et al. (2014) and Schaye 
et al. (2014) does not discern between radio and quasar 
mode. A further difference between Illustris and Booth & 
Schaye (2009) is that they use a different accretion model. 
Instead of a constant value for the accretion enhancement, 
a is set to 1 at low gas densities and varies as a power law 
of density at higher densities (see Booth & Schaye 2009, for 
a discussion). In addition to this, Schaye et al. (2014) in¬ 
clude a model which can lower the accretion rate in order 
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Figure 9. Spatial region corresponding to the ‘ejected baryons’ 
range of Table 2. In addition to the dark matter density cut, we 
here require a temperature higher than T > 6 x lO'^K. 



Mpc 
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Figure 10. Temperature field of the gas in the slice of Fig. 3. 
We see that the gas forms hot bubbles around the haloes and 
filaments. 



Mpc 


Figure 11. Plot of the gas metallicity for the slice which has been 
used in Fig. 3. The extended baryon distributions which protrude 
into the dark matter voids are rich in metals, suggesting that 
feedback processes are the cause of the extended distribution. 


to account for the angular momentum of accreting gas (see 
Rosas-Guevara et al. 2015). Also, they use a radiation effi¬ 
ciency of cr = 0.1 and a coupling efficiency of Cf = 0.15. In 
order to prevent that the injected energy is cooled away im¬ 
mediately, the feedback energy is not emitted continuously, 
but only once enough energy has accumulated to lead to a 
heating AT of 10® K. The energy in then injected into one 
particle close to the black hole (or multiple particles if the 
energy produced per time-step is high enough to heat several 
particles to 10® K). 

There are also differences with respect to the used hy- 
drodynamical scheme. Vogelsberger et al. (2012) found that 
the SPH code GADGETS produces higher temperatures in 
haloes than the moving-mesh code AREPO, which has been 
employed in Illustris. They argue that this higher tempera¬ 
ture has to do with spurious dissipation in SPH related to 
the higher numerical noise in this method. This might partly 
explain why the AGN feedback efficiency parameters had to 
be set higher in Illustris than in a comparable SPH simula¬ 
tion in order to reach the same suppression of star formation 
in massive haloes. Also, the mass resolution in McGarthy 
et al. (2010) and Le Brun et al. (2014) is very low compared 
to Illustris, which could lead to different density estimates 
at the black hole and thus influence the accretion rates sig¬ 
nificantly, even if the same accretion model was used. 

For lower mass haloes AGN feedback is less important, 
as lower mass haloes typically have less massive black holes 
and the AGN luminosity scales with the square of black 
hole mass. Therefore, at low to medium masses the dom¬ 
inant feedback process regulating star formation is stellar 
feedback. A part of the stellar feedback is already implic¬ 
itly included through the effective equation of state for star 
forming gas (Springel & Hernquist 2003). In addition to 
this, Illustris uses a non-local ‘wind’ feedback scheme (see 
Vogelsberger et al. (2013)), where a part of the star form¬ 
ing gas can be stochastically converted into wind particles. 
These wind particles are allowed to travel freely for some 
time (or until a density threshold is reached) before their 
mass and energy is deposited into the appropriate gas cell. 
Without these feedback processes, galaxies would form too 
many stars (the so called over-ooling problem). Thus the 
parameters of this models will also influence the results we 
presented (especially the region below 10^^ M 0 in Fig. 1). 
The Illustris stellar mass function is a bit too high in this 
mass range, suggesting that the feedback should be more ef¬ 
fective in preventing star formation. Similarly, we find a high 
baryon fraction around unity in these haloes. As stated in 
Section 3, there are substantial observational uncertainties 
for the baryon fraction in this mass range, and the obser¬ 
vations we included in Fig. 1 are to be thought of as lower 
limits. 

We should caution however that while the feedback of 
Illustris is so strong that it depletes the massive haloes of 
their gas, the stellar mass function of Illustris is still too high 
for haloes with a stellar mass above 10^^M© (see Fig. 3 in 
Genel et al. 2014). This indicates that the feedback is not 
efficient enough in suppressing star formation at this mass 
scale. Simply tuning the energy parameters of the feedback 
model cannot resolve the tension between the low gas con¬ 
tent in massive haloes on one hand and the too high stellar 
mass function on the other hand. Either the AGN feedback 
model has to be modified (Genel et al. 2014, suggests making 
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Figure 12. Redshift evolution of the contribution of haloes, voids and filaments to the total baryonic mass (a) and the total dark matter 
mass (b). While the evolution of the baryons and the dark matter is similar at redshifts higher than two, at low redshifts the haloes, lose 
baryons due to feedback processes while the dark matter in haloes is still growing. 



Figure 13. Evolution of metal fractions in gas and stars. The 
values are normalized to the total metals at 2 : = 0, and the con¬ 
tribution of haloes, filaments and voids to the metals budget in 
gas is shown. 


the duty cycle longer and thus avoiding short strong bursts) 
or this is a hint that other, so far neglected processes are 
important. It will be interesting to see if improved feed¬ 
back models, which reproduce the baryon fraction in mas¬ 
sive haloes correctly, will continue to find such a significant 
amount of baryons in dark matter voids. 


5 SUMMARY 

We investigated the global mass distribution of baryons and 
dark matter in different structural components of the Illus- 
tris simulation. By attributing dark matter density regions 
to the constituents of the cosmic web, we could measure the 
mass and volume of haloes, filaments and voids, as summa¬ 
rized in Table 2. 

At redshift 2 : = 0 we find that the haloes host 44.9 % 
of the total mass, 49.2 % of the dark matter and 23.2 % of 
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Figure 14. Mass fraction of gas to total baryons in temperature 
versus baryon density space. We subdivide the histogram into 
hot gas, diffuse gas, condensed gas and warm-hot intergalactic 
medium (WHIM). The mass fractions in these regions at 2 : = 0 
are indicated in the figure. 


the baryons. Their volume fraction amounts to only 0.16 % 
of the total simulation volume. The filaments host 44.5 % 
of the dark matter and 46.4 % of the baryons and corre¬ 
spond to 21.6 % of the simulation volume. The dark matter 
voids contain 6.4 % of the dark matter but 30.4 % of the 
baryons, and the volume fraction of the voids amounts to 
78.2 %. About 80 % of the baryons inside voids correspond 
to gas which has been ejected from haloes due to feedback 
processes. Most of this material, which occupies one third of 
the volume of the voids, has been ejected at redshifts lower 
than z = 1. 

An analysis of a comparison run, in which star forma¬ 
tion, cooling and feedback had been switched off, showed 
that the baryons trace dark matter well in a non-radiative 
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Figure 15. Redshift evolution of the mass fraction in the gas 
phases of Fig. 14. In addition to the gas phases, we also include 
the contribution of stars. 

model, which is not the case for the full physics simula¬ 
tion. This demonstrates that the feedback model has signif¬ 
icant consequences for the global mass distribution. While 
the feedback model of Illustris succeeds in reproducing a 
realistic stellar mass function, the low baryon fraction in 
massive haloes suggests that the feedback model is not cor¬ 
rectly reproducing their gas content. At the moment it is still 
unclear whether our analysis would produce similar results 
with a different feedback model. Work on improved feed¬ 
back models is ongoing, and it will be interesting to clarify 
this question in future studies based on a next generation of 
simulation models. 
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